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We study the electronic properties of an array of small metallic grains connected by tunnel 
junctions. Such an array serves as a model for a granular metal. Previous theoretical studies of 
junction arrays were based on models of quantum dissipation which did not take into account the 
diffusive motion of electrons within the grains. We demonstrate that these models break down at 
sufficiently low temperatures: for a correct description of the screening properties of a granular metal 
at low energies the diffusive nature of the electronic motion within the grains is crucial. We present 
both a diagrammatic and a functional integral approach to analyse the properties of junction arrays. 
In particular, a new effective action is obtained which enables us to describe the array at arbitrary 
temperature. In the low temperature limit, our theory yields the correct, dynamically screened 
Coulomb interaction of a normal metal, whereas at high temperatures the standard description in 
terms of quantum dissipation is recovered. 

PACS numbers: 73.40Gk, 72.15.-v, 74.50+r 



I. INTRODUCTION 



During the past two-decades much attention has been devoted to the study of disordered bulk metals with electmru 
electron interactionsErO from the one hand and arrays of normal metallic grains connected by tunnel junctionsETtl 
from the other. The Coulomb interaction plays a very important role in both types of system and many non-trivial 
physical effects occur due to it. In disordered metals, e.g., the Coulomb interaction results in quantum corrections 
to conductivity and reduces the density of states, ajJieaomenon known as zero bias anomaly. Qualitatively similar 
effects were predicted^ and studied experimentallvQoaEIl for arrays of tunnel junctions in the normal state, where a 
sufficiently strong Coulomb interaction can suppress the conductivity (Coulomb blockade of tunneling). 

In spite of the similarities in the physical behavior of disordered metals and normal arrays of tunnel junctions, 
completely different approaches are used for the theoretical description of these respectiiie. systems. As concerns 
disordered metals with interaction, either diagrammatic methodsEJ or non-linear CT-modelsEiFt£l are employed. The 
starting point of these theories is a model of a weakly disordered bulk metal. On the other hand, for the description 
of Coulomb blockade effects in a junction array o^a granular material, a completely different approach based on the 
theory of Ambegaokar, Eckern and Schon (AES)EJ is commonly used. The free energy functional proposed by AES 
does not contain any disorder: the finite conductance obtained within this model is due to tunneling of electrons 
between grains and can be incorporated into the theory through a term describing the so-called quantum dissipation, 
first introduced phenomenologically in Ref. [l^. 

Although both approaches have been introduced almost 20 years ago, only a few-indirect attempts have been made 
to reconcile them. We mention the work by NazarovEil and Levitov and ShytovO that treaJ_.thc diffusive zero bias 
anomaly and the Coulomb blockade of tunneling on the same footing. More recently, NazarovE3 showed how Coulomb 
blockade phenomena survive in diffusive systems. Indeed, the properties of a granular metal without electron-electron 
interactions are qualitatively similar to those of a disordered bulk metal (see e.g. Ref. |l9| ) and it would be natural to 
expect that this qualitative similarity persists even in the presence of interactions. In fact, one may conjecture that 
bulk disordered metals from the one hand, and granular metals or arrays of tunnel junctions from the other, can be 
described within one unifying scheme. 

In order to develop such a scheme, we consider in this paper a granular normal metal with Coulomb interaction 
between electrons. The consideration is simplified by taking into account only the long range part of the interaction 
that leads to Coulomb blockade effects. We assume that the macroscopic dimensionless conductance qt of the 
granular metal is large, which enables us to develop a perturbative theory with g^, 1 as the small parameter. Within 
this theory, physical quantities of interest can be calculated at arbitrary temperatures. The diagrammatic analysis 
is supplemented by the derivation of a cr-model that can be considerably simplified provided one does not take into 
account weak localization effects. 
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We demonstrate below that granular metals and networks of tunnel junctions can always be described using the 
standard techniques developed in the theory of disordered metals (expansions in cooperons and diffusons, a- model 
calculations). With the help of these techniques one can in principle calculate any physical quantity without using the 
notion of quantum dissipation. Disorder is inevitably present in the system even if electrons move ballistically within 
the grains. In the latter case any small irregularity in the shape of the grains causes the electron motion within the 
grains to be chaotic and this assumption is sufficient for our theory to be applicable. An integrable shape of the grains 
does not seem to be realistic, and even if it happened for isolated grains, tunneling from grain to grain would addjaa 
additional chaoticity. These ideas about intrinsic disorder in granular metals have been used in previous studiesEjEil. 

We find that, at not too low a temperature, the AES free energy functional gives correct results. However in the 
limit T — > AES theory is no longer applicable. The region of validity of the AES free energy functional is determined 
by the inequality T g^S, where 5 = (vV) 1 is the mean level spacing, v is the density of states at the Fermi surface 
and V is the volume of a single grain. 

We conclude in particular that the AES notion of quantum dissipation at T = can only be applied to tunnel 
junctions connecting infinitely large conductors. If the volume of the conductors is finite the AES picture can be 
applied at finite temperatures only. Apparently, the notion of quantum dissipation should be treated with care. A 
qualitative discussion of the problem "dissipation versus dynamical screening" in the case of superconducting grains 
can be found in a recent lecture by AltshulerE3. 

The remainder of the paper is organized as follows. In Sec. |l| we formulate the model and write physical quantities 
in terms of functional integrals over anticommuting fields. We decouple the Coulomb interaction by an integration 
over auxiliary fields. In Sec. Ill we consider the AES action and show that it does not correspond to the physics of 
normal metals with a screened Coulomb interaction. In Sec. IV we develop a diagrammatic technique. We show what 
kind of diagrams correspond to the AES action and demonstrate that, at low temperatures, additional contributions 
arise. Summing up these additional diagrams we arrive at expressions corresponding to the dynamically screened 
Coulomb interaction of a normal metal. Sec. ^ contains a derivation of a new action for a granular metal with 
electron-electron interaction which is applicable at arbitrary temperatures. The results are compared with the results 
obtained diagrammatically and with those obtained on the basis of the AES action. Our results are summarized in 
the Conclusion. 



II. CHOICE OF THE MODEL 



We consider an array of normal metallic grains coupled to each other. The dimensionality of the array may be 
arbitrary. The grains are assumed to contain imperfections: there can be impurities inside them as well as on their 
surface. This implies that the electron motion in the grains is chaotic. The mean level spacing in a single grain is 
equal to 6 = (fV) 1 . However, the shapes of the grains need not considerably differ from each other and we assume 
for simplicity that the grains are arranged in a cubic lattice. The electrons can hop from grain to grain; moreover, 
they interact with each other. We are interested only in the long range part of the Coulomb interaction and write it 
in a simplified form describing charging of the grains. In such a formulation the spin of the electrons is not important 
and can be taken into account at the end of the calculations when writing proper densities. 

Under these assumptions the electron motion can be conveniently described by a functional integral of the type 



exp(-SM)Xty, 
with an action S [ip] that can be written in the form 



The action S s [tp] in Eq. (2.1) describes the electron motion within the grains and we write it in the form 

SM =Y,J Q drdT ^ ^ r ) ( d r-V + Hi) Mr, r), (2.2) 

where ipi (r, r) is a fermion field in grain i at imaginary time r and coordinate r (tp* is its complex conjugate), /i is 
the chemical potential, and (3 — T" 1 is the inverse temperature. The fermion fields ipi (r) satisfy the condition 

ipi (t) = -i>i (t + 0). (2.3) 

The operator Hi in Eq. (|2.2|) includes scattering by impurities and can be written as 
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Zm 



(2.4) 



where m (r) is the impurity potential. We assume that it is Gaussian distributed with the correlation 

1 



(ui(r)uj(r')) 



S(r - r')6 i: 



imp 



The action St [ip] stands for tunneling between the grains, 

S tbP]=y2 d Ti } i(T~)tijll>j(T), 
ij J ° 

and S c [ip] describes the charging of the grains 

ij ° 

where 

m (t) = / drip* (r, r)%pi (r, r) 



(2.5) 



(2.6) 



(2.7) 



is the density field in grain i, e is the electron charge, and Cij is the capacitance matrix. Eqs. (2.1 - 2.7) describe the 
model completely and one can start explicit calculations. 

Due to the Coulomb interaction, the action S[ip] is not quadratic in ip and the integration over ip cannot be 
performed immediately. A convenient way to proceed in such a case is to decouple the term S c [ip] by a Gaussian 
integration over auxiliary fields. This transformation can be written as follows: 



V(iP\iP)e~ SW =N VVe~ S2[v] / V{iP* ,iP)e~ Sl[,l) > v] 



(2.8) 



where TV is a normalization factor and Vi (r) is the decoupling field. The action S2 [V] in Eq. (2.8) has the form 



(2.9) 



which shows that the variable Vi (r) has the meaning of a voltage on grain i at the time r. The new effective action 
Si [ip, V) reads 



S[tP,V] = S g [iP,V] + S t [iP], 



with 



S g [iP, V\=Y.[ dT ^*( r ) { d - + + ieVi{r)) ^(r). 



(2.10) 



One more transformation can be performed exactly. Following the procedure of Ref. where a single grain was 
considered, we represent Vi (r) in the form 



Vi (r) - y 4 ° + V (r) , 

where V t ° is the static part and V, (t) dr = 0. The replacement 



ipi (r) -» ipi (t) exp {-i(j>i (r)) 



(2-11) 



(2-12) 



with 



(r) = e / V (r') dr' 



(2-13) 
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does not violate the condition given by Eq. (2.3) and we can remove the variable Vi (r) from the action S R [if),V], 
Eq. ( ]2.10| ). However, this variable appears in the action St[ip] as an additional phase of the field tp (r). As concerns 
the static part V®, its fluctuations can be neglected even in a single isolated grain, provided the temperature is high 
T <5cj. This restriction is not necessary for the system of the coupled grains in the limit of large conductance 
gr 3> 1 conside red here. 

[2.13 ) we finally reduce the calculation of physical quantities to the computation of a functional 



Using Eqs. Q2.1 



integral of the form 



exp (-SbM - S^j, 4>] ~ S a [<f>]) £>W, 



(2.14) 



where 



So VIA =Y, ( T ) ( 9 r " /* + H i ) ' > < " ) ' 



(2.15) 



Sify, 4>]=^2 dTtijil>*(T)il>j(T) exp (icj) tj (r)) . 



(2.16) 



Jo 



Cy #i(r) #j(r) 



2e 2 



(2.17) 



and <^>ij ( t) = <fe (t ) - ^ (r). 

Eqs. (|2.14 ■ 2.17 ) completely specify the model that will be studied in the subsequent sections. Disorder is still 
present in the operators Hi, Eq. (2.4), but the mean free path within the grains is assumed to be large, Iko 3> 1, 
where fco is the Fermi momentum. In this limit, the macroscopic conductivity is determined mainly by the tunneling 
conductance gr, 



g T = ir/(2e 2 R T ) = 2tt 2 *V « 6A5kn/R T 
where Rt is the tunneling resistance. 



(2.18) 



III. QUANTUM DISSIPATION DESCRIPTION OF GRANULAR METALS WITH COULOMB 

INTERACTION 

A functional integral formulation was used in Ref. ^ to treat the quantum dynamics of a Josephsorujunction. The 
action derived in that work (AES action) was used later in the context of a normal tunnel junctional, as well as in 



connection with arrays of tunnel junctions in a number of paperstl. The AES action can be obtained from Eqs. (2.14- 



2.17). First, one should average the Green function Go corresponding to the action Sq[iJ)], Eq. (2.15), over impurities 



in the first Born approximation, which gives (in Fourier representation) 



<V,,(pi=li e -tfp) + i^M] . CU) 



2Ti m p 



where T lmp is defined i n Eq . (2.5) and £(p) = p 2 /2rn — /i. After that one should make a cumulant expansion in 
the term Si[ip,<p], Eq. ( |2.16| ). Keeping only the second order of the expansion and assuming for simplicity that the 
capacitance matrix Cy is diagonal one arrives at the AES action S*aes ot— i 



/ \ 2 

Saes^H^E / dT (^) +Y E / rfrrfr'a(r -r')(l- cos (^(r)-^(r'))), 
1 o iJ o 



(3.2) 



where 



a{T-r') = -^- -. (3.3) 

sm irl (r — r) 
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In Eq. (3.2) E c is the charging energy E c — e 2 /2Cu. The second term in Eq. (3.2) contains the sum over neighboring 
grains and t is th e tunneling energy between them. In the limit of small phase fluctuations we can expand the second 
term in Eq. (3.2) with respect to faj and obtain 



S[<f>] 



1 o IJ o 



)(fa j (r)-fa j (r')T 



Using the periodicity of all functions in /3 we write the Fourier expansion as 

ipij (t) = T}^ (iu; n ) exp(-iw„r), a(r) = T 2J a(icu n ) exp(-iw Tl r), 



(3.4) 



(3.5) 



where lu u — 2imT are Matsubara frequencies. Hence we obtain for the action 

sy>] = ^E T E w « i i 2 +?E T E i ^ ii M*»n) i 2 



(3.6) 



The second term in Eq. ( |3.6| ) describing the tunneling is linear in frequency and keeps its form down to T = 0. 
Therefore, it was attributed to quantum dissipations. Performing a Fourier transformation in space 



i{uj n ) = ^0 Wn (k)exp(ikRi), 



(3.7) 



we rewrite the action S\d] as 



S[4>] 



^E 




^(l-coskd Q ) |^„(k)| 



(3.8) 



where d a , a 
function 



x, y, z are vectors connecting the centers of the neighboring grains. The phase-phase correlation 



n 



AES 



K,k) =T(^ B (k)0* (k)> 



can be immediately calculated using Eq. (3.S) and we obtain 



IIaes (w n ,k) = 



4B„ 



f ffT £; c |w„|^ Q (l-coskd a ) 



(3.9) 



Suppose that one assumes, following Refs. that Eq. ( |3.2| ) is applicable to a system of tunnel junctions down to 
T = 0. This would imply in particular that Eq. (3.9) is correct for T — > 0, too. However, the derivative d^i/dr is 
proportional to the voltage on the grains which in turn is linearly related to the charge density. Multiplying Eq. ( |3.9| ) 
by u^, one thus obtains the density-density correlation function K (w„,k). Let us consider limit of small frequencies 
and wave vectors k. In this limit, the network should correspond to a bulk normal metal, where on e sh ould have 
the propagator corresponding to a screened Coulomb interaction. However, the limit uj, k — » in Eq. (3.9) does not 
reproduce the screened Coulomb interaction. 

How to resolve this discrepancy? In the next sect ion we carry out an expansion in and demonstrate that 
keeping only the second order, which led us to Eq. (3.2) is not sufficient and that, in the limit of small frequencies or 
temperatures, one should sum up an infinite class of diagrams. In section [v] we will then show how these processes 
can be included into an effective action formulation structurally similar to the AES approach. 



IV. SCREENED COULOMB INTERACTION IN GRANULAR METALS 



In this section we present a detailed derivation of the density-density correlation fcmction K (w n ,k) summing an 
infinite series of diagrams. For the computation we use the standard diagram techniquccil with a modification suggested 
in Refs. |2^,^l] to include tunneling between the grains. Within this technique one should make an expansion in the 
tunneling t{j, the phases fa and in the impurity potential Ui(r). Below we denote the electron Green functions 
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by solid lines, the pha ses <pi by wavy lines and the tunneling elements by crossed circles. Impurity propagators 
corresponding to Eq. (^4) are denoted by dashed lines. 

Before turning to the quantitative discussion, let us briefly outline the skeleton of the analysis. To lowest order 
in the tunneling matrix elements, the polarization operator of the system is represented by the process depicted in 
Fig. |^. This diagram describes the (absolute square of) amplitude for tunneling from one grain into the next and the 
subsequent relaxation inside the 'target grain' (represented through the Green function) . The quantitative evaluation 
of this process leads to the coupling constant gr of the AES approach. 

On this level, coherence effects associated to mulitiple tunneling and/or impurity scattering are completely neglected; 
the tunneled electron 'forgets' about its phase memory implying that the dissipative picture of the AES approach 
obtains. Indeed, the two Green functions depicted in Fig. || live in different grains which means that no phase 
correlation is possible. How can this picture change in principle? Including one more order in the tunneling, leads to 
corrections of the structure shown in Fig. ||. What makes these processes qualitatively different from those discussed 
above is that they include the possibility of phase coherent multiple impurity scattering: The two Green functions 
labeled i or j scatter off the same impurities and/or chaotic potential where the associated scattering phases may 
cancel due to the fact that one of the Green functions is advanced, the other retarded. The net two-particle modes 
emanating from such processes, 'diffusons' in the context of disordered systems, are represented through the hatched 
regions of Fig. |^, right. The key question now is, are such processes relevant or not? A crude estimate can be given 
as follows: It is known that for low temperatures (for a precise definition of 'low', see below), each diffuson mode 
diverges as S/T. The fact that we are considering processes of higher order in the tunneling amplitude introduces one 
more power of <?t- Thus, for T <~ 5qt, higher order processes become as important as the first order contribution. 
For lower temperatures, these processes must be taken into account to obtain a physically correct picture. In the 
following we will put this discussion onto a firm basis and discuss how the coherent tunneling series can be summed 
in a controlled way. 

Before starting the computation for a granular system we want to mention an important difference between the 
diagrams for the current-current and density-density correlation functions. Calculating the classical conductivity for 
a granular system in the lowest order in the tunneling we need not renormalize the current vertices by impurities. 
To understand this fact let us recall a well known result for bulk metals. The diagram for the Drude conductivity is 
shown in Fig. [|. 





FIG. 1. Diagram for the Drude conductivity and current vertex for a bulk metal. 



The second diagram in Fig. renormalizing the current vertex, gives an additional integral over momentum that 
yields zero, 



- /pGo £ (p)G 0e (p)d 3 p = 0. 
m J 



(4.1) 



In Eq. (4J), the Green function Gq £ (p) is determined by Eq. (3^). Diagrams with many impurity lines yield zero 
as well and the current vertices are not renormalized. The same is true for granular metals, which can be seen after 
a proper replacement of the current 



eid sinpd. 



(4.2) 



As concerns the polarization operator or the density-density correlation function, the situation is different. In this 
case, instead of vector vertices we have scalar ones and the renormalization due to impurities can be important. 

Expanding the functional integral, Eq. ( 2T4 ) in the tunneling term Si, Eq. (2.16), we obtain to second order in iy 
and in (pi the diagrams represented in Fig. |2| 




FIG. 2. Polarization operator for a granular metal in the lowest order in tunneling. 



G 



The two electron lines relate to different grains. As the impurities in the different grains are not correlated, impurity 
lines connecting two Green functions are do not exist (although the diagram represents a polarization loop.) The 
second and third diagrams in Fig. || are equal to each other and have an opposite sign with respect to the first diagram. 
The analytical result for the sum of the three diagrams reads 



P (w„,k) = -\wn\(2/ir)g T ^2(l - cos k a d) , 



(4.3) 



where the dimensionless tunneling conductance gr is given by Eq. ( 2.18 ). We see from Eq. (4.3) that, although 
the Green functions Go E (p) determined by Eq. (3.1) contain Ti mp , all information about the disorder drops out of 
the polarization loop. To calcu late the phase-phase correlation function one should notice that the bare propagator 
corresponding to 5*2, Eq. (2.17), has the form 



n (w„,k) = — 5- 

and diverges at to — > 0. Therefore, one should sum up the infinite geometrical scries represented in Fig. 



(4.4) 





FIG. 3. Diagrammatic analogy of AES action in the limit of small phase fluctuations. 



The summation can be easily performed and we arrive at the propagator Haes (w,k), Eq. (3.9). Summarizing , we 



conclude that, in order to reproduce the results obtained in the quadratic approximation of the AES action, Eqs. (3.4, 
3.G, 3.S), one should sum up the diagrams of Fig. ^. On the other hand, the form of the propagator IIaes ( w n> k) does 
not correspond to the propagator of a normal metal with a screened Coulomb interaction (see e.g. [l]), which signals 
that something is missing in this analysis. 

To resolve this paradox one should consider diagrams of higher order in tunneling, which complicates the calcu- 
lations. First of all, let us note that in order to consider higher order diagrams one should specify tunneling more 
explicitly. We assume that the radius of the area at which the grains contact each other much exceeds atomic distances 
(large number of conducting channels). At the same time, the potential barrier between the grains can be arbitrarily 
large, which can lead to an arbitrary tunneling conductance. Such a situation can be conveniently modeled by random 
tunneling matrix elements t pq correlated as 



(tpqtp' qf) — t (Spp' Sqq' 



&pq' &p' q) 



(4.5) 



where t pq are written for the same contact 



Correlations of the matrix elements of different contacts are assumed to 
be zero. To simplify notations we do not write subscripts numerating the contacts. Taking into account Eq. (4.5) 
is indeed very important. In order to illustrate this fact we consider two diagrams represented in Fig. |Ja. We see 
immediately the difference between them: the first diagram is not averaged over t pq and can contains four different 
momenta. After averaging over t pq one obtains the second diagram containing only three different momenta. 

As it has been mentioned, we will expand not only in the tunneling amplitude t but also in the phases 4>. Performing 
the expansion in the phases and denoting them by wavy lines we have three different possibilities as depicted in Fig. 
|b. 

a) 





b) 



Pi 



HP 



V \s — ^- 



'l a 



vP 



vPi 



1-4 r 4 14 

FIG. 4. All classes of diagrams which contribute to the polarization operator of granular metals in the fourth order of 
tunneling. 
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Now we can start averaging over impurities. As usuafia, we average first the Green functions and obtain Eq. (3.1) 
for the average Green function G$ s (p). Then, we can consider averages with impurity lines connecting the Green 
functions on opposite sides of the diagrams in Fig. |i|b. Diagrams with non-intersecting impurity lines are most 
important and we recognize the well-known diffusons. 

In the present paper we consider the limit of sufficiently small grains, such that all relevant energies like temperature 
T, frequency u>, tunneling energy t, etc. are much smaller than the Thouless energy Et = 7r 2 Dq / R 2 , where Dq is the 
diffusion coefficient of a single grain and R is the radius of the grain. In this limit the spectrum of the diffuson in 
a single grain is zero-dimensional (no dependence on momenta inside the grain). Coupling between the grains leads 
to a dependence on quasi-momentum k. Diagrams that should be summed in order to give the complete form of the 
diffuson are represented in Fig. 0a. 



a) 




imp 




complete 
diffuson 



b) 




c) 



1 J 
— ® — 

— <S>— 

1 j 



1 J 1 



+ 



+ 



— ® — <S> — 
i j i 



FIG. 5. Equation for complete diffuson of granular metals. 

Now let us present some details of the calculation of the diagrams in Fig. [|. The different classes of diagrams 
contributing to the polarization operator of the granular metals before impurity averaging are shown in Fig. ^b. 
Before impurity averaging the sum over Matsubara frequencies for the first diagram in Fig. Wb is 



T^G(ie n + iw„,pi)G(ie„ + iui n , p 2 )G(ie n , pi)G(ie n , p 4 ). 



(4.6) 



Writing this expression we used Eq. (4.5). After impurity averaging of the first diagram in Fig. ^b, we obtain the 
diagram represented in Fig. ||a. 

Now the disorder averaging can be done and we obtain, in a standard way, the average of the product of two Green 
functions for a granular metal, (G R G A ) — D(u) n ,q). Here 



D = +D^ZD, 



(4.7) 



(uj n , k) — 2nv/\uj n \ being the diffuson f or a single isolated grain. The equation for the complete diffusonEll is 
shown in Fig. ||b. The self-energy £ in Eq. (4.7) is the sum of the three diagrams in Fig. |5|c. The second and third 
diag ram s in Fig. ||c are equal to each other and have the opposite sign with respect to the first diagram. Solving 
Eq. (F7) we obtain 



2-kv 



+ {2/n)g T 5 £ (1 - cos k a d) 



(4.8) 



The calculation of the second diagram in Fig. |]b is analogous to the calculation of the first one. The third diagram 
in Fig. ^b is independent ofthe frequency u> n and yields a constant. The diagram in Fig. ^|a contains only one diffuson. 
More complicated diagrams can be drawn and one of them is depicted in Fig. ^. 
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FIG. 6. Example of a high order correction to the polarization of granular metals. 



However, all such many-diffuson contributions can be neglected provided the tunneling conductance <?t, Eq. (2.18), 
is large . A direct calculation of this diagram shows that it involves the small parameter 1/<?t- This is similar to 
expansions in the diffusion modes for a bulk metal where the inverse conductivity is the expansion parameter. The 
condition g? ^> 1 means that we are far from the Anderson metal-insulator transition. 

The result of averaging and summation of the diagrams represented in Fig. can be written as 



Pi(w n ,k) 



cu n \5 ( (2/7t)3tI](1 - cosk a d) 



\<jj n \ + (2/ir)g T 8 XX 1 ~ cos k a d) 



(4.9) 



We note here that the functions Po and Pi hav e a different sign. Comparing the function Pi (oj„,k), Eq. (4.9), with 
the quantum di ssipa tion part Po (<^ri,k), Eq. (4J3), represented in Fig. || we see immediately that the contribution 
Pi (cj„,k), Eq. (4.9), can be neglected only if the temperature T is sufficiently high T <y \w n \ 3> gr$ (we recall that 
the static component of the phase <f> has been neglected and therefore u> n ^ in Eq. (4.9). So, the notion of quantum 
dissipation is applicable essentially at non-zero temperatures only for a granular material. However, it can be used 
for a contact connecting two bulk metals with an infinite volume. 

At low temperatures T < gx8, there is no reason to neglect the contribution Pi (w„, k), Eq. ( [4.9[ ), and its presence 
changes completely the form of the propagator. Adding the functions Po (u> n , k), Eq. (fi~3|), and Pi (u> n , k), Eq. (p~s|), 
we obtain the total polarization P (u) n , k) 



\uj n \ 2 (2/ir)g T J2{l - cosk a d) 



P(Wn,k) 



(2/ir)g T 8 XX 1 ~ cos k ad) ' 

a 

Proper diagrams contributing to this function are represented in Fig. ^. 



(4.10) 




J 



+ ( in M{ + (all permutations) 

FIG. 7. Relevant polarization loops for a granular metal. 



imp. 



Now, we can write the final expression for the phase correlation function n {u> n , k) that should be written instead of 
IIaes (^n: k), Eq. (|3.9|). Using the bare propagator n (w„, k), Eq. (4.4), and the self-energy part P (oj n , k), Eq. (4.10), 
we write the phase propagator II (u/„, k) in the form 



n(w„,k) = 



n„ 



4P C /^ 



i -n n p 



(2/7T)g T J2( 1 - cos k * d ) 
1 + AE C | (i , ri | + ( 2 / ff ) gr 5 ^(i_ cos fe a <i) 



(4.11) 



Using the relation between the phase 
Coulomb interaction 



<p and the voltage V, Eq. (2.13), we derive an equation for the effective 
V eH (u> n ,k) =«4n(w„,k). (4.12) 



On the other hand, the dynamically screened Coulomb interaction in a disordered bulk metal can be written in the 
form 
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wher e Vp is the bare Coulomb potential. We see that Eqs. ( 4.11 . 4.12 ) written in the continuum limit correspond to 
Eq. |Ll~3l ). 

Thus, we conclude that, when calculating physical quantities, the main contribution is due to small k only in the 
limit T <C <?t<5- Provided this inequality is satisfied, the properties of the tunnel junction array are equivalent to 
the properties of a weakly disordered bulk metal with Coulomb interaction. In the opposite limit T 3> grS, one can 
describe the system in the "quantum dissipation" approximation. 

Calculations were performed in this section using a diagrammatic expansion in the lowest order in the field c6. On 
the other hand, the AES free energy functional, Eq. ( |3.2| ), is written for arbitrary values of (f>. Can one generalize 
Eq. ( 1.12] ) to arbitrary values of <p1 We will try to do this in the next section, deriving a simplified version of a replica 
a- model and using a non-trivial saddle point. 



V. EFFECTIVE ACTION FOR GRANULAR METALS WITH COULOMB INTERACTION. 



In the preceding sections we analyzed effects of the Coulomb interaction on the electron motion in granular materials 
using a diagrammatic expansion. This approach works well when fluctuations of the phase <f> or, in other words, 
fluctuations of the voltages on the grains are small. However, one can imagine situatipas when the fluctuations are 
large. For example, one can try to consider non-perturbative excita tion s like instantonsE3 corresponding to an integer 
charge transfer between grains. In principle, the AES action, Eq. (3.2), can describe such excitations very well but, 
as we have seen previously, it can only be used at sufficiently high temperatures T 3> gT$- So, our task now is to 
derive an action that will be applicable at lower temperatures. 

The derivation at temperatures T < gxS is more complicated because now we should explicitly take into account 
disorder, which was not necessary for the derivation of the AES action (strictly speaking, disorder determines the 
mean free time Ti mp in the Green functions used for derivation of the AES functional but this is a trivial contribution). 
On the other hand, the problem is not as complicated as the problem of Anderson localization in the presence of 
interactions considered in Refs. In our diagrammatic expansions we considered the limit of large conductances 

gT 3> 1 and neglected all weak localization corrections. 

Nevertheless, we startirfche derivation of the action using the replica cr-model approach for interacting systems 
suggested by FinkelsteinErliil. Necessary simplifications will be made in the cr-model. Although there are difficulties in 
using the replica cr-models for non-perturbative calculations, our goal is more modest and the final results obtained 
in the limit gx 3> 1 neglecting all weak localization corrections will not depend on replica indices at all. 

The model under co nsiderat ion has been formulated in Section ||. The derivation of the cr-model can be carried 
out starting from Eqs. (2.1-2.7). Instead of the fields ip and V one should write fields ip — {ip a } and {V a } carrying 
the replica index a = 1, ....r. Upon calculating physical quantities for an arbitrary r, one should put r — 0. 

A proper cr-model for bulk systems with interaction has been derived in Refs. |3|JTl|. A generalization to granular 
metals can been made without difficulties following Ref. |l^, where a cr-model was written for a granular metal 
without interaction. Although one can proceed along the lines of Refs. i.e. expanding the action obtained after 

integration over ip in the field V, a more economic way is to use the substitution given by Eqs. (2.12, 2.13), as well as 
the subsequent Eqs. (2.14-2.17). Of course, this is already an approximation because we neglect the static component 
of the phase <p. However, for our purposes this is not an essential restriction because we do not want to consider 
effects of localization. Moreover, working in the limi t T 6 o r gr 3> 1 allows us to ignore the static component of 6. 

The derivation of the cr-model starting from Eqs. (2.14-2.17) is practically the same as the one presented in Ref. [l|. 
As usual, one decouples the "effective interaction" of the type ip A that appears after averaging over impurities by a 
Gaussian integration over matrices Q. These matrices contain as elements both th e rep lica indices and Matsubara 
frequencies (or two imaginary times r and t'). Neglecting the tunneling term, Eq. (2.16), one would obtain a zero- 
dimensional cr-model (we assume that all relevant energies are smaller than the Thouless energy Et of a single grain) . 
The relevant tunneling term in the cr-model is obtained by a cumulant expansion in the tunnelin g. Th e only difference 
with respect to the model without interaction of Ref. [l^ is the presence of the phases (j> in Eq. (2.16), which leads to 
additional phase factors in the term describing the coupling between the grains. 

The final result for the free energy F[Q, <fi] can be written as 



F[Q,( 



Fi 



F^+Ft[Q,> 



(5.1) 
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dfr (t) dtj (r) 
dr dr 



dr. 



F*lQ] = jY,[ tr(<D T Qi(r,T))dT, 

F T [Q,0]=-^X; I tr (e^WQi (r,r') e^^Qj (r',r)) drdr', 



(5.2) 
(5.3) 
(5.4) 



where (^>y (r) = ^ (r) — (f>j (r) and the symbol tr implies a trace over both replica indices. The field <fi is diagonal in 
the replica indices, while the matrix Q (t, t ) is a 2r x 2r matrix with the constraints 



J Q (r, t") Q (t", t') dr" = 8 (r - r') , trQ (r, r) 



= 0. 

The action of the operator cD in Eq. ( |5.3| ) on an arbitrary function / (r, r) is given by the relation 

d 



uf (r, t') 



-i lim --/ (r, r') 



(5.5) 



(5.6) 



In frequency representation this ope rator is equal to a vector of fermionic Matsubara frequencies. The matrix Q (t, t'), 
satisfying the constraints, Eq. (5J)), can be conveniently parametrized as 



(r,r')= / U (t,t") A t „ iT ,„U (t"' , /) drV , 



where 



A T , T / = - 



iT 



sinT7r (r — r') 



and J7, U are unitary matrices 

/" U (r, t") [7 (t", t') dr" = S (r - t') . 
The matrix A in frequency representation has the form 

Amn — ^mn^gn (Cn) • 



(5.7) 



(5.8) 



(5.9) 



The limiting cases of the free energy F[Q, cj>], Eqs. (5.1-5.6) are rather simple. Without electron-electron interaction 
the phase difference between neighboring grains equals to zero 4>ij = 0. Such an action for granular metals corresponds 
to that derived within the supersymmetry scheme by one of the authorsEj. The second limiting case is achieved at 
sufficiently high temperat ures . In this limit disorder is not import ant and the matrix Q does not fluctuate, being 
equa l to the value A, Eq. (5JJ). Inserting Q = A into Eqs. (5^1-^4) we reproduce immediately the AES action, Eq. 
( |3.2| ) (in fact we obtain the r times repli cated A ES functional but only one replica field is important for us). 

At not too high a temperature, Eqs. (5.7-5~9|) lead to a non-trivial behavior even if the disorder is weak and all 
effects related to weak localization can be neglected. The limit of weak disorder enables us to neglect fluctuations 
of Q at given <fi and take its value from a saddle-point equation. We thus reduce the computation of the functional 
integral with the free energy functional to the solution of a saddle-point equation, substituting the solution to the free 
energy F[Q, </>], Eqs. (5.1-5.4), calculating a functional integral with the reduced free energy F[4>] 

F[4>] = F[Q, 4>] (5.10) 

where Q is the solution of the saddle-point equatio n for a given ip. 

Minimizing the functional F[Q, </>], Eqs. (5.1-5.4), with the constraints, Eqs. (5.5), we obtain the following equation 

I W^ (T) Q l (r,r")e-M-"),Q j (r",r')}dr" + %[Q% (t, t) ,u>)] = 0, (5.11) 



9r 
2 



where the summation over j is performed over the nearest neighbors of i . In the commutator [. . , . .] one should exchange 
properly the times when changing the order of the functions. 
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The Eq. (5.11) represents an integral equation that enables us to find, in principle, the solution Ki (r) for arbitrary 
values of the parameter gTS/T. In general, this equation can be solved only numerically. However, for gx 3> 1 drastic 
simplifications arise. For large intergranular coupling, the phase 4> fluctuates only weakly from grain to grain which 
means that fcj is small. Under these conditions, a solution of the diagonal form 



Q, 



iKiir') 



(5.12) 



can be sought for. As the solution is assumed to have a diagonal form, all replica indices decouple and t he equ ations 
can be solved for each replica separately. Therefore we can drop the replica indices and consider Eq. (5.11) as an 
equation for a single function Ki (r) for a given function <pij (r) . 



In the high temperature limit T 3> grS the second term in Eq. (5.11) is much larger than the first one and we 
arrive at the solution 



Ki (r) = 0, 



(5.13) 

which leads to the AES free energy functional, Eq. ( |3.2| ). 

In the opposite limit of low temperatures T <c grS Eq. (5.11) can be simplified using the assumption that Ki (r) 
varies slowly in space. Then, one may expand the exponentials in Kij (r) = Ki (t) — Kj (t), which leads to the 
equation 



g T MiKj + <fe) ) A] + j[^-,K i ] = 0. 



(5.14) 



Eq. ( 5.14 ) is presented without writing the integration over imaginary times r explicitly, but of course this integration 
is implied. Transforming this equation to frequency and momentum space we obtain 

^2g T ^(l -cosg/d)(l - sgn(cj n .)sgn(w m )) + ^(w„ - w m )(sgn (u n ) -sgn(u> m )J K q ^ rn = (5.15) 



-—V q , n - m (sgn (ui n ) - sgn (u> m )), 



which gives the solution 



Kk.n = —Vk. n i- 



sgnw„ 



(2/ir)g T 8 ]C (1 - coskd a ) 



(5.16) 



where Kk,n and Vk,n — (i/ e)u) n <j>nk are the Fourier transforms of the functions Ki (r) and Vi (r) respectively. Substi- 
tuting this ansatz for K back into the action we obtain 



F[V] 



1 



l + 4E r 



(2/7r)gTEa( 1 - cosfc a d) 
| + (2/7r) ffT( 5^ Q (l-cosfc a d) 



(5.17) 



where E c is the charging energy. In the limit of small characteristic momentum, k <C we recover the diffusively 
screened Coulomb interaction in a disordered evironment. 

From here one can proceed to the calculation of physical observables such as the conductance, the density of 
states or others. To this end one should introduce configurations Q = T QT fluctuating around the mean field 
configuration discussed above. Next one would follow the standard algorithm of doing calculations within the a- 
model approach: (i) express the quantity of interest in terms of Q, (ii) compute the functional average over the 
effective action F[Q] — F[V,T] as good as is possible. 



Let us outline how the connection between the enlarged AES-type formulation, (5.1) - (5.4) and the Finkelstein 
approach for interacting disordered media can be made explicit. To do so, we subject our Q's to a gauge transformation, 



J4>, 



-i<j>i 



(5.18) 



as a result (i) the hopping part of the action, Ft becomes (^-independent . However (ii) in the frequency part, 
uj t — > o) T + d T <f>i. Gaussian integration over 9^ then produces an effective action F[Q] which, after taking a continuum 
limit, is identified as the action of the Finkelstein approach. 
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VI. CONCLUSION 



In this paper we constructed a theoretical framework to describe electronic transport in arrays of tunnel junctions 
or granular metals, in the limit of large inter-granule conductance, <?t- Both disorder and the electron-electron 
interaction were taken into account. The prime objective of this enterprise was to unify two large, and seemingly 
non-overlapping theories of interacting metallic compounds: the AES approach, focusing on charging phenomena 
in individual grains and their impact on large scale transport behaviour, and the Finkelstein theory of disordered 
interacting metals with its emphasis on the interplay interaction/disorder. The key to the reconciliation of these 
two approaches was to observe that for low enough temperatures, T < gxS, the effective action underlying the AES 
theory becomes incomplete. The physical reason is that for low temperatures the electrons and holes participating 
in the tunneling processes between individual grains maintain their quantum phase memory for a long time, largely 
in excess of the average tunneling time. This means that the particle tunneling is not only accompanied by charging 
and dissipation (as in the AES approach), but also by other, more long ranged physical processes. Specifically, we 
identified the quantum dissipation contained in the AES action as_the high temperature limit of the long ranged 
screened electron-electron interaction, an observation first made inE3. We re-emphasize that these phenomena were 
not bound to the presence of a significant disorder concentration; indeed, none of the results discussed above, displayed 
dependence on some 'disorder concentration'. The only thing that mattered was chaoticity of the electron motion on 
scales set by the phase coherence length, a condition that is practically always met in real life systems. 

Technically, two different routes for including these processes into the theory were proposed. One consisted of the 
perturbative summation of diagram classes associated to multiple chaotic scattering. The other, based on an effective 
action formulation, indentified what is 'missing' in the AES action. Remarkably, this second approach readily led 
to a unification of the AES approach and the Finkelstein nonlinear er-model for weakly disorderd interacting metals. 
With the benefit of hindsight, this fact is easy to understand: It is a well known fact in mesoscopic physics, that the 
long processes resulting from multiple scattering can be interpreted as a certain type of Goldstone modes. Within 
an effective action approach, these modes must be described by an own degre e of f reed om, the Q-matrices of the 



nonlinear cr-model. Thus it is no surprise that we obtained an action of the type (5.1 - p.4| ), involving Q-fields and the 
Coulomb phase fields of the AES approach, as an effective description of the low energy phase of the system. For large 
temperatures, the fluctuations of the Q-matrices became inessential and the AES action was retrieved. In contrast, 
the continuum limit of this action was identified as the Finkelstein model and the Q-matrices smoothly integrated 
into the AES approach, without leading to conceptual complications. 

Let us make some remarks on potential experimental ramifications of our findings. It has been predicted in 
Refs. ^]|7] that two-dimensional normal arrays of tunnel junctions should undergo a Kosterlitz-Thouless-Berezinskii 
(KTB) phase transition at a temperature T c of the order of the charging energy E c . At low temperatures, T < T c , the 
array should be in an insulating state (Coulomb blockade), whereas for T > T c the array is conducting. In particular, 
at T > T c , the conductivity a of the array is predicted to increase with temperature according to a square-root cusp 
dependence, cr(T) ~ exp(— 2b/ y/T/T c — 1), where b is a constant of order unity. The experimentsLjntrEiJ mentioned 
in the Introduction searched for this KTB transition. In view of our results, the applicability of the KTB scenario 
should be governed by the parameter g^S/T,,. The experiments of Refs. [7| ]9|]l0| w ere done on arrays with relatively 
large grains, such that gxS/T c <C 1. Indeed, the results obtained in Refs. |7|,[10| were in agreement with the KTB 
scenario. However, Ref. || found thermally activated behavior of cr(T), rather tharuthe predicted square root cusp 
dependence on temperature. On the other hand, the experiment by Yamada et al.u was done on arrays consisting 
of relatively small Cu grains with a size of about AO A, such that gxS/T c ~ 20. Clearly this is beyond the range of 
applicability of Refs. |j7|and the KTB scenario-should be treated with care in this case. Remarkably, Ref. || concludes 
good agreement of the results with the theorytnl Given these discrepancies, we conclude that a careful analysis of the 
temperature-dependent conductivity of junction arrays in the framework of the theory presented in this paper would 
be of interest. 
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